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O ' In "extended phase space" approach to quantum geometrodynamics numerical solutions 

to Schrodinger equation corresponding to various choice of gauge conditions are obtained 
for the simplest isotropic model. The "extended phase space" approach belongs to those 
Q^i appeared in the last decade in which, as a result of fixing a reference frame, the Wheeler 

! - DeWitt static picture of the world is replaced by evolutionary quantum geometrody- 

namics. Some aspects of this approach were discussed at two previous PIRT meetings. 
We are interested in the part of the wave function depending on physical degrees of 
freedom. Three gauge conditions having a clear physical meaning are considered. They 
are the conformal time gauge, the gauge producing the appearance of A-term in the 
^ I Einstein equations, and the one covering the two previous cases as asymptotic limits. 

tJ" ' The interpretation and discussion of the obtained solutions is given. 
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1. Introduction 

In this paper we present solutions to quantum geometrodynamical Schrodinger equation cor- 
responding to various choice of gauge conditions for the simplest isotropic model. It is widely 
accepted in quantum geometrodynamics to illustrate general ideas taking simple cosmological 
models as examples. The reason why physicists working in this field appeal to simple models 
is that now quantum geometrodynamics is just as far from being a completed theory as it was 
decades ago. One must confess that hitherto there is no agreement on what "first principles" 
this theory should be based and what is the form of master equation for a wave function of the 
Universe. The first version of quantum geometrodynamics, proposed by Wheeler and DeWitt 
[HE], encountered a number of fundamental problems (for discussion, see |3l|4l|5]). The main 
problem is the so called "frozen formalism", or the absence of time evolution. It is easy to 
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see that the source of the problem of time consists in the apphcation of the Dirac postulates 
to gravitational field, according to which not the Schrodinger equation but the constraints as 
conditions on a wave function play the central part in the theory. As a result of impossibility 
to resolve the problems of the Wheeler - DeWitt quantum geometrodynamics in its own limits, 
in the last decade there appear a new tendency in the development of the theory which can 
be called Evolutionary Quantum Gravity. This tendency may be characterized by the two fea- 
tures: firstly, the recognition of the fact that it is impossible to obtain the evolutionary picture 
of the Universe without fixing a reference frame and, secondly, the rejection of the Wheeler - 
DeWitt equation and the reestablishment of the role which the Schrodinger equation plays in 
any quantum theory. 

The tendency embraces several approaches (see, for example, P,[7j, where a dust fluid is 
considered as a good choice to fix a reference frame in quantum gravity), to which the "extended 
phase space" approach belongs. Some aspects of the latter were discussed at two previous PIRT 
meetings [HI E]. The approach is based on a careful analysis of peculiarities of quantization of 
the Universe as a whole [TOl [11]. The analysis showed that quantum geometrodynamics as a 
mathematically consistent theory failed to be constructed in a gauge invariant way, therefore, 
the Wheeler - DeWitt equation, being a constraint on a state vector, loses its significance and 
should be replaced by a gauge dependent Schrodinger equation resulting from the Hamiltonian 
formulation of the theory in extended phase space. A wave function satisfying the Schrodinger 
equation is determined on extended configurational space that involves gauge gravitational 
degrees of freedom equally as physical ones. However, we are actually interested in the part 
of the wave function depending on physical degrees of freedom only, since this very function 
defines probability distributions of physical quantities. 

In Section 2 we shall describe the model and the Schrodinger equation for the physical part 
of wave function for the given model. Since the form of the Schrodinger equation is gauge 
dependent, to obtain descriptions of the Universe corresponding to various gauge conditions 
(in other words, to various reference frames) one has to solve, in fact, absolutely different 
differential equations. It naturally leads us to the question, is there any correspondence among 
solutions of the equations? And how should they be interpreted? 

Let us note that while in [6l [7] the authors work with a certain parametrization of gravita- 
tional variables (as a rule, it is the Arnowitt - Deser - Misner parametrization [121) some 
"privileged" reference frame, our approach, though was applied to cosmological models with 
finite degrees of freedom, aimed at including arbitrary parametrizations and a wide enough 
class of gauge conditions. We shall consider three gauge conditions having a clear physical 
meaning: the conformal time gauge, the gauge producing the appearance of A-term in the Ein- 
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stein equations, and the one covering the two previous cases as asymptotic hmits. For a closed 
universe, the first and third gauges gives rise to a discrete Hamiltonian spectrum, while the 
second gauge leads to a continuous spectrum. From a pure methodical viewpoint, the first and 
third cases are much easier to be treated, and in Section 3 numerical solution for these cases 
will be presented, meantime the second case admits qualitative consideration only. Section 4 
contains physical interpretation and conclusions. 



2. The model and the Schrodinger equation for the physical part 
of the wave function 



The action for a closed isotropic universe is 

' 1 ad? 
2~W 2 



S^^at) = -Jdt Na'e{a) , S^gf^ = J dt ttq (^N - ^c^j. (2.2) 

Matter fields are described in this model phenomenologically, without a clear indication on the 

nature of the fields. The dependence of its energy density e{a) on the scale factor a determines 

its equation of state, namely, for the power dependence e(a) = — , the equation of state is 

a" 

known to be P{mat) = ("3 " ^{mat), £0 is a constant whose dimensionality in the Plank units 
is ppilpi- Since we are interested in early enough stages of the Universe evolution, we shall 
suppose that the Universe was filled with radiation with the equation of state p(mat) = i^^imat), 
i.e. 

e{a) = ^. (2.3) 

S(gf) is a gauge- fixing part of the action, its variation giving rise to gauge dependent terms in the 
Einstein equations. In ordinary quantum theory this terms are to be excluded by asymptotic 
boundary conditions. As was argued in [10], in the case of the Universe with a non-trivial 
topology, which, in general, does not possess asymptotic states, making use of asymptotic 
boundary condition is not justified. 

If so, the gauge-fixing action describes a subsystem of the Universe, some medium, whose 
state is determined by a chosen gauge. In (12.21) a differential form of the gauge condition 

N-f{a)=0 (2.4) 

is used. The equation of state for this subsystem is 

1 f'{a) 

P{obs) = - j.^^^ a£(^obs)- (2.5) 
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The index (obs) indicates that this subsystem corresponds to an observer studying the Universe 
evolution in his reference frame. 

The action (12.11) is a particular case of the action for a cosmological model with a finite 
number degrees of freedom considered in [8|, [11] . The Schrodinger equation for the physical 
part of the wave function looks like 



--\ r \ r- + -Na^ - iVa^£(a)^ 

2\ a da \^ a da j 2 

From the classical point of view, E is given by 



= E^!. (2.6) 

N=f{a) 



E = -jV^Tl^^,,,)d'x. (2.7) 

'^fl{obs) is ^ quasi energy-momentum tensor obtained by variation of the gauge-fixing action; it 
is not a real tensor in the sense that it depends on a gauge condition. T^(oft^) describes the 
subsystem of the observer in the gauged Einstein equations [8]. It can be shown that the 
integral (12.71) of T^f^obs) taken over space is a conserved quantity for the class of gauge conditions 
(12. 4p . Thus, E characterizes the energy of the observer subsystem. 

It may be said that on a phenomenological level this approach takes into account interaction 
between the observer subsystem and the physical Universe. The interaction causes rebuilding of 
energy balance of two subsystems. It is expected that at the late stage of the Universe evolution, 
when the Universe is well described by General Relativity, gauge effects are negligible, and the 
values of E must be very close, if not equal, to zero. However, at the early quantum stage E 
may have essentially non-zero values, and the exploration of its spectrum is the main task of 
this work. 

Now we consider several gauge conditions. 

1. The conformal time gauge N = a. The equation of state of the observer subsystem is 
the same as that of the matter: p[obs) = ^^^(obs)- Substituting N = a and (12.31) in (12. 6p . we get 

E + eo^E (2.9) 

lrf2^ 1 



After redefinition 



we obtain the equation 



+ -a^^^J = E^. (2.10) 
2 rfa2 2 ^ ^ 

Therefore, Eq. ( ]2.10p describes the Universe filled with a "substance" with the equation of 
state p = -6. Just some part of the energy of this substance may be due to a usual matter 
while the other part may be due to gauge, or observer, effects. 
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It was shown in [13] that Eq. fl2.10l) can be obtained in the limits of the Wheeler - DeWitt 
quantum geometrodynamics by rewriting of the Wheeler - DeWitt equation H"^ = as a 
Schrodinger-like equation H'^ = E^^. Under additional requirements, that imply choosing a 
certain gauge condition and including a certain kind of matter into the model, the classical 
Hamiltonian constraint H = can be presented in a new form, H = E, H = H — E, where 
E is a. conserved quantity which appears from phenomenological consideration of this kind of 
matter. So, in this approach, E = Eq, i.e. E is entirely due to the usual matter (radiation). 

On the other side, the need for making a choice of gauge to rewrite the Wheeler - DeWitt 
equation in the special form H"^ = E'^ witnesses to gauge noninvariance of the Wheeler - 
DeWitt theory. As was already emphasized above, the Wheeler - DeWitt equation loses its 
meaning, and it seems to be reasonable rejecting it rather trying to hold it by any means. 

The effective potential U{a) = -c? is given at Fig. 1(a). 
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Fig.l. The effective potentials for Eqs. (|2J^ . (|2JJ]) . (f2A3l) . 

2. Na?" = 1. The gauge is believed to produce the appearance of A-term in the Einstein 
equations, since it is the analog of a more general condition det||5f'^'^|| = 1. The equation of 
state p(obs) = —£(obs)- The Schrodinger equation takes the form 

7—^ + —— + —*-^^ = ^^. (2.11) 

2 da? a^ da 2a^ a^ ^ ' 

Here Eq characterizes a contribution of the matter fields (radiation). If one includes into the 

model de Sitter false vacuum with the equation of state P(o6s) = —£(ohs) and the dependence 

E{a) = Eo, it does not affect the form of the equation ( 12. lip after redefinition ( 12. 9p . Then one 

could say that vacuum energy as well as gauge effects are responsible for eigenvalues of E. 

The effective potential U(a) = — - — ^ depends on the parameter ^o- According to modern 

2a^ a'^ 

cosmological notions, the Universe was created in a metastable state under the barrier depicted 
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at Fig. 1(b) and then tunneled through the barrier. The smaller the parameter Eq is, the higher 
and narrower the barrier becomes. There is a non-zero probability for arbitrary large values 
of the scale factor a; it means that the Universe may expand to infinity in spite of the sign 
"+" we have put before the second term in (12. ip . which corresponds to the closed model. It 
demonstrates that a naive correspondence between the kind of a cosmological model and the 
form of the effective potential has no grounds. 

3. = a-\ — -. This gauge covers the two previous cases as asymptotic limits. The equation 
of state is 

1 a'' — 3 m Tn\ 

P{obs) = - ^4 ^ ^ ^(ofes) • (2.12) 

At a — ^ the equation gives P(obs) = —^(obs)', at a ^ cxd it gives P(obs) = 7t^(o6s)- Again, after 
redefinition (12. 9p the Schrodinger equation looks like following 

2 V aV da^ da 2 2a^ ^ ' 

It is easy to check that Eqs. (12. lip . (]2.10p are the asymptotic limits of (I2.13P at a ^ 

and a ^ oo respectively. In this case the Universe is believed to be filled by some mixture 

of matter and vacuum. In consequence of the redefinition (12.90 . the value of E is due to 

matter contribution as well as gauge effects. Like in a previous case, the effective potential 

Uia) = -c? H — ^ depends on the parameter eo and depicted at Fig. 1(c). The barrier 

2 2a^ 

at small a disappears when £0 = and Eq > 0.1. The potential for some value of Eq is shown 
at Fig. 2. One can see that the potentials of Eq. (12. lip (green graph) and of Eq. (I2.10p (blue 
graph) are asymptotic forms of the potential of Eq. ( I2.13p . 
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Fig. 2. The effective potentials for some values of Eq. 
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3. Numerical solutions 



The Hamiltonian operators in Eqs. fl2.10l) . f l2.13p have a discrete spectrum, and one meets no 
technical difficulties to obtain numerical solutions to these equations. The operator in fl2.6p is 
Hermitian for an arbitrary gauge (12.41) if the measure in Hilbert space of solution is taken to be 

M(„) = (3^1) 

One can see that the measure, like the equation itself, is gauge-dependent. 

The standard method of finding eigenvalues and eigenfunctions consists in the expansion 
onto a basis functions which are orthonormal on the interval [0, oo] with the measure (13. ip : 

vI/(a) = Ec„^:(a); (3.2) 

n 

= \/7^S^=7T^'^^^") = \iu^^ i^-^YaiL^M: (3.3) 
V in + s)\ ^M{a) \ {n + s)\\ a J 

oo 

J ij'^*{a) M{a) da = 6nm, (3.4) 



L'^{a) are Laguerre polynomials. The problem is reduced to finding eigenvalues and eigenvectors 
of the Hamiltonian matrix in the basis (13.30 . The more terms are held in the expansion (13. 2p . 
the higher the precision is. The results of calculations of ffist five eigenvalues are presented at 
Table 1. 



Table 1. 



Eq. (12.101). TV = 


a 


1.5 


3.5 


5.5 


7.50001 


9.50008 






= 


2.87886 


5.32668 


7.66977 


9.9591 


12.2175 




eo = 


1/500 


2.87846 


5.32635 


7.66947 


9.95882 


12.2173 




eo = 


1/150 


2.87754 


5.32558 


7.66877 


9.95817 


12.2166 




= 


1/50 


2.87489 


5.32337 


7.66677 


9.9563 


12.2149 


Eq. rt2.13ll. 


Bo = 


= 1/2 


2.77519 


5.24152 


7.59315 


9.88783 


12.1496 


N = a + \ 


eo 


= 1 


2.66102 


5.15088 


7.51266 


9.81349 


12.0792 




So 


= 3 


2.04887 


4.72486 


7.14847 


9.48369 


11.7714 




So 


= 4 


1.59368 


4.47069 


6.94071 


9.29951 


11.602 




Bo 


= 5 


0.972188 


4.1924 


6.71849 


9.1044 


11.4236 




Bo 


= 7 


-1.07592 


3.59902 


6.25063 


8.69468 


11.0497 
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One can see that for Eq. (12.101) . N = a the spectrum is equidistant, the difference between 
eigenvalues is equal to 2 in the Plank units (the deviation from this value is entirely due to 
calculation inaccuracy). 

In the case of Eq. (I2.13p . N = a-\ — -, the eigenvalues do not differ significantly for Eq < — 

a-^ ^ 50 

and converge to limiting values at = 0. For eq > — the spectrum levels tend to go down 

50 

into the potential pit. The schematic picture of the spectrum is shown at Fig. 3. 
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Fig. 3. The spectrum levels for some potentials. 



Fig. 4-6 pictures the probability distributions for the first (ground state), third and fifth 

solutions to Eq. (I2.10p and Eq. (I2.13P when = and Eq = 7. One can see that at 

150 

the qualitative level the probability distributions do not significantly differ. The peak of the 
probability distribution in the all cases tends to shift to large values of the scale factor a 
for larger eigenvalues of E. One could expect this result since the matter and gauge effects 
contribute to the value of E. So, when the energy of matter increases, there may be enough 
probability for the scale factor to reach large values. 
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4. Concluding remarks 

We should recognize that we have considered a very simple model and the obtained results are 
not of high degree of generality. The present work is just a small step "to find the way". 

We have seen that the second gauge condition, = — , leads to a continuous spectrum 
of eigenvalues of the Schrodinger equation (12. 6p . while the two other gauges, N = a and 

= a H — -, leads to a discrete spectrum, in other words, the second case is substantially 
different. It seems that one should seek for the reason in the structure of spacetime. Indeed, 
the gauge N = — corresponds to the Universe in which the interval of proper time between 
two subsequent spacelike hypersurfaces tends to zero as a — oo, meantime it is not the case for 
the two other gauges. Since any gauge condition determines the form of the effective potential, 
this circumstance require a more careful exploration. It would be interesting to study the gauge 
A^ = 1 H — r, for which at a — oo the reference frame becomes a synchronous one (A^ = 1) and 
the equation of state of the observer subsystem at a — > oo is that of dust: P(obs) = 0. 

The resemblance of probability distributions for solutions to Eqs. (12.101) . (12.131) also deserves 
our attention. It demonstrates that one can reveal some relation among solutions for certain 
classes of gauge conditions. Let us note that this problem is well-known in the Wheeler - 
DeWitt quantum geometrodynamics, and the question how solutions to the Wheeler - DeWitt 
equation are related, was discussed as soon as parametrization noninvariance of this theory had 
been realized. Then Halliwell [H] proposed to restrict the class of admissible parametrizations. 
Since parametrization and gauge conditions have a unified interpretation [5], it implies also a 
restriction of the class of admissible gauge conditions, i.e. such an approach implies that it is 
permissible to describe the Universe in only one or several "privileged" reference frames. This 
way seems to be artificial since we do not know for sure what reference frame is privileged. Our 
point of view is that we face a new problem of finding classes of gauge conditions within which 
solutions to the Schrodinger equation are stable enough with respect to a choice of gauges, 
the determination of the classes seems to be inseparable from our understanding of spacetime 
structure. 
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